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ABSTRACT 

Context. TW Hya is a classical T Tauri star that shows significant radial-velocity variations in the optical regime. These variations 
have been attributed to a 10 M Jup planet orbiting the star at 0.04 AU. 

Aims. The aim of this letter is to confirm the presence of the giant planet around TW Hya by (i) testing whether the observed RV 
variations can be caused by stellar spots and (ii) analyzing new optical and infrared data to detect the signal of the planet companion. 
Methods. We fitted the RV variations of TW Hya using a cool spot model. In addition, we obtained new high-resolution optical & 
infrared spectra, together with optical photometry of TW Hya and compared them with previous data. 

Results. Our model shows that a cold spot covering 7% of the stellar surface and located at a latitude of 54° can reproduce the 
reported RV variations. The model also predicts a bisector semi-amplitude variation < 10 m/s, which is less than the errors of the RV 
measurements discussed in Setiawan et al. (2008). The analysis of our new optical RV data, with typical errors of 10 m/s, shows a 
larger RV amplitude that varies depending on the correlation mask used. A slight correlation between the RV variation and the bisector 
is also observed although not at a very significant level. The infrared H-band RV curve is almost flat, showing a small variation (<35 
m/s) that is not consistent with the published optical orbit. All these results support the spot scenario rather than the presence of a hot 
Jupiter. Finally, the photometric data shows a 20% (peak to peak) variability, which is much larger than the 4% variation expected for 
the modeled cool spot. The fact that the optical data are correlated with the surface of the cross-correlation function points towards 
hot spots as being responsible for the photometric variability. 

Conclusions. We conclude that the best explanation for the RV signal observed in TW Hya is the presence of a cool stellar spot and 
not an orbiting hot Jupiter. 
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1. Introduction 

' TW Hya (spectral type K7, M, =O.7M , age~10Myr) is 
an extensively studied cl assical T Tauri star (CTTS) (e.g. 
Rucinski & Krautter|[l98ll) . With a distance of 56+7 pc (ESA 
19971) . it is probably one of the closest laboratories for studying 
of our own solar system in its early stages of formation. 

Recently ISetiawan etaTJ ([2008. SHL08, hereafter) an- 
nounced the discovery of a giant planet (~10Mj up ) orbiting TW 
Hya at a separation of 0.04 AU. Their discovery is based on high- 
resolution optical spectroscopic observations: the analysis of two 
independent datasets obtained with a difference of ~ 3 months 
shows a significant radial velocity (RV, hereafter) variation with 
an amplitude of 196+61 m/s and a periodicity of 3.56~0.02 days. 
This, together with the lack of correlation between the RV vari- 
ation and the cross correlation-function bisector (BIS), has been 
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interpreted as proof of the existence of a planetary-mass object 
orbiting the star. 

In this letter, we study whether the RV variations in TW Hya 
can be caused by stellar spots and not by a planetary-mass com- 
panion. With this in mind, we have tried to reproduce SHL08 
data with a cool spot model. In addition, we obtained almost si- 
multaneous optical & near-infrared high resolution spectroscopy 
and optical photometry of the source in order to detect the planet 
signal with a different dataset. 

2. A spot model for TW Hya 

To test whether a spotted surface is a viable model for explain- 
ing the SHL08 dataset, we simulated the changes of TW Hya's 
line profile due to such a spot. For that purpose, we used Soap 
(Bonfils & Santos 2008, in prep.), a program that calculates 
the photometric, radial-velocity, and line-shape modulations in- 
duced by one (or more) cool stellar spot. 

Soap computes the rotational broadening of a spectral line by 
sampling the stellar disk on a grid. For each grid cell, a Gaussian 
function represents the typical line of the emergent spectra. It has 
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Fig. 1. Radial velocity and bisector variations as calculated by 
our spot model. In the top panel, the two epoch (clear and dark 
symbols) RV measurements from Setiawan et al. (2008) are su- 
perimposed. 

the same width, <x , as the typical spectral line of a star with no 
rotational velocity (or one too low to be resolved) and of similar 
spectral type than the star of interest. The Gaussian is Doppler- 
shifted according to the projected rotational velocity (v sin z) and 
weighted by a linear limb-darkening law (a = 0.6). For a given 
spot (defined by its latitude, longitude, and size), Soap computes 
which of the grid cells are obscured and removes their contribu- 
tion to the (integrated) stellar spectrum. Finally, the stellar spec- 
trum output by Soap is the sum of all contributions from all grid 
cells. In our case, the stellar spectrum is an averaged spectral 
line. It behaves like the cross-correlation profile of stellar spec- 
tra, and Soap uses it to compute the RV and the bisector, in ad- 
dition to the photometric variation. 

In the case of TW Hya, we chose /?* = 0.9 R (follow- 
ing SHL08) and an inclination of i = 7°, whi ch is the likel y 
inclination of the TW Hya circumstellar disk dOi et alj 12004). 
The strong veiling in the atmosphere of TW Hya prevents us 
from using its B-V color index as a temperature proxy. Both cr 
and v sin i will probably yield wrong estimates if derived from 
a calibration based on this color index. Therefore, we only used 
the B-V color to derive a first estimation of cr = 4.45 kms" 1 
dMelo et al.ll200ll) . We then iterated on different values to show 
that our conclusions are unchanged. The rotational period of the 
star, the spot size and latitude, and the star velocity (y) were left 
as free parameters. Adjusting these parameters, we were able to 
reproduce SHL08's RVs with a dark spot covering about 7% of 
the stellar surface, located at a latitude of 54° and with a rotation 
period of 3.56 days. With the assumed stellar radius and inclina- 
tion, and with the fitted rotational period, we derived a v.sini of 
1.56 kms" 1 . 

Remarkably, our solution has x 2 - 3.27, very similar to the 
Keplerian solution derived by SHL08 (^ 2 =3.32). Our model also 
indicates that such a spot would produce a 4% photometric vari- 
ability in the V-band, and a peak-to-peak amplitude change <10 
ms~' in the BIS measurements (see Fig. 1). An extreme choice 
of cro =1 kms -1 increases the amplitude of the BIS variation up 
to <30 ms _1 . 

SHL08 measured the BIS and conclude that the observed 
RV change is not due to a spot. However, their BIS measure- 
ments have individual errors of ~50-100 m/s and their overall 



dispersion is almost as large as the RV variation, therefore much 
higher than the variation expected from our model. Varying the 
constraint on the inclination, we find solutions with^ 2 <3.6 from 
i = 2° to i ~ 80°. With cr =4.55 kms -1 , the peak-to-peak varia- 
tions of BIS were found to be lower than 100 m s _1 for all incli- 
nations below ~ 50°. 

According to our model, the low inclination of TW Hya's 
system weakens significantly the power of a bisector analysis 
and other diagnostics should be used to validate or disprove the 
planet hypothesis. 



3. New observations of TW Hya 

3.1. High-resolution optical spectroscopy 

With the goal of confirming the period, phase, and amplitude of 
the RV signal observed by SHL08, we used the CORALIE spec- 
trograph (1.2-m Euler Swiss telescope, La Silla, Chile) to ob- 
tain a series of 12 consecutive radial-velocity measurements of 
TW Hya. The data were gathered between February and March 
2008. Each measurement has an accuracy of about 10 ms -1 . 
CORALIE has a spectral resolution R-A/AA ~ 50 000, and is 
known to have a long term precision in radial- veloc ity of a few 
~3ms-' (e.g. lMavor et~aill2004tlWilson et al.ll2008l) . 

The CORALIE RVs are obtained by cross-correlating the 
stellar sp ectrum with an optimized weighted template (for de- 
tails see iPepe et aTll2002l) . For TW Hya, we performed correla- 
tions using cross-correlation masks that are appropriate for stars 
with different spectral types: G2, K0, K5 (the one closer to the 
spectral type of TW Hya), and M4. This allows to verify that 
the amplitude of any putative RV signal shows the same value, 
independently of the mask used. This is expected if the radial- 
velocity variations observed stem from the presence of a planet. 
On the other hand, an amplitude dependence on the mask used 
could be expecte d if the RV signa l is produced by the presence 
of a stellar blend (Santos et al. 20Q2|) or a cool stellar spot. In the 
latter case, the correlation with a mask built for a later spectral 
type (e.g. M4) should provide a signal with a smaller amplitude 
than the one derived using a mask for an earlier spectral type 
(e.g. G2). This is due to the fact that later type masks are more 
sensitive to the lines in the spot than are the earlier type masks. 
In the former case, the spectroscopic contrast between the spot 
and the rest of the photosphere is thus weaker. 

The CORALIE RVs have been complemented by older RVs 
measurements of TW Hya (three datapoint s) obtained with the 
HAR PS spectrograph (data described in Guenther & Es positol 
120071) . and the FEROS RVs listed by SHL08. Using the HARPS, 
FEROS, and CORALIE (K5-mask) radial-velocities together, 
we find a clear signal with a period of ~3. 56 days, and a semi- 
amplitude of K=238+9ms~ 1 (Fig. [2] upper-right panel). The 
parameters of a Keplerian fit are presented in TableQ] Using 
only the points listed by SHL08 and our own fitting rou- 
tines, we do confirm their solution, with a semi-amplitude of 
K=196+13 m s _1 . If caused by the presence of a planet, our data 
could thus imply a planetary mass that is 20% higher than the 
one announced by SHL08. 

Interestingly, however, a look at the CORALIE measure- 
ments alone shows that by using different cross-correlation 
masks we can derive significantly different amplitudes when 
fitting a Keplerian function to the data (Fig. [2] ; bottom-right 
panel). For instance, using a mask built for a G2-type dwarf 
we derive a value of K=292+35ms~ 1 . This value decreases to 
K=274+29ms- 1 , K=252+20ms- 1 , and K=169±21 ms" 1 , us- 
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Fig. 2. Top left: Phase-folded radial-velocity measurements of 
TWHya with the Keplerian fit obtained with the HARPS 
(squares), CORALIE (filled circles), and FEROS (open circles) 
data. The best Keplerian fit to the three groups of points is 
shown. The semi-amplitude of this fit is K=238 m s _1 . Top right: 
Original data from SHL08. The best Keplerian fit is shown 
(K=196ms _1 ). Mid and lower left panels: CORALIE phase- 
folded radial-velocities derived using three different CCF masks 
(see text for details). The same Keplerian function as found us- 
ing SHL08 is shown for comparison. Lower right: RV amplitude 
of the best Keplerian fit to the CORALIE data derived using 4 
different masks. 



ing the KO, K5, and M4 masks, respective^. This is in line with 
the presence of a spot instead of a planet to explain the observed 
RV variations. 

If the signal is caused by a stellar spot, a color effect can also 
be expected. If we correlate the spectra with a given mask using 
only the blue or the red sides, we can expect to obtain slightly 
different values for the RV amplitude. This difference comes 
from the different contrasts (flux ratio) coming from the spot 
and the surrounding photosphere as a function of wavelength. In 
a simple blackbody approximation, the contrast between a spot 
and the surrounding disk is stronger at bluer wavelengths, thus 
implying a stronger RV signal. 

To test this, we derived CORALIE RVs using the K5 mask 
but correlating only a red portion of the spectrum (4150-4850A) 
and a blue portion of the spectrum (6000-6700A). In the for- 
mer case, we obtained a semi-amplitude of K=228+25ms~ 1 , 
while in the latter we derived K=210+20m s _1 . These values are 
marginally compatible with a small variation. Indeed, a simple 
comparison of the effect due to the different flux ratios in the dif- 
ferent bands used (considering a spot with a temperature around 
3000 K and the TW Hya temperature of 4000 K - SHL08) shows 
that a variation of about 10% can be expected. Following this ra- 
tionale, the RV amplitude in the infrared (IR) regime is expected 



1 To derive these values we fixed the eccentricity to zero. No major 
differences were found if the eccentricity is set as a free parameter. 



200 


_ . 1 i i i i 1 < i i i 1 i i < i 1 i . 






1 » 


;^\+ + 






s 

-100 






- 


-200 


Slope=-0.3±0.1 
, 1 1 , , 




+: 




-200 -100 100 

VRAD [m/s] 


200 





Fig. 3. BIS vs. radial-velocity (derived with the K5 mask) for the 
CORALIE data of TW Hya. For clarity, the vertical and horizon- 
tal scales were set to be the same. The best fit to the data is also 
shown. 



to be much smaller (~5 times smaller) if the RV signal is indeed 
caused by a spot. We test this hypothesis in Sect. 3.2. 

Table 1. Elements of the fitted Keplerian function. 
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As explained in Sect. 2, the absence of any clear correla- 
tion between the RV variation and the bisector of the cross- 
correlation function was considered by the SHL08 as evidence 
of a planet orbiting TW Hya. However, the errors in their data 
(around 100ms 1 ) may have hidden any low amplitude signal. 
To test this, we computed the BIS v alues for our CORA LIE 
measurements (for the definition see lOueloz et al.l 1200 ll) . As 
can be seen in Fig. [3] the BIS and RV are slightly correlated 
(slope=-0. 3+0.1), though not at a very significant level (cor- 
relation coefficient, r=0.82). This result is compatible with the 
output from the simple model described in Sect.|2] 

Finally, we searched for phase variations in the RV signal be- 
tween the FEROS (SHL08) and CORALIE data, obtained with a 
time span of about one year. The results show that no measurable 
phase shift was found over more than 300 days to a precision 
ofl%. 

3.2. High resolution near-IR spectroscopy 

A planet orbiting TW Hya should be detectable with a similar 
RV amplitude and orbital period at all wavelength ranges. As 
mentioned above, if the RV signal reported by SHL08 is caused 
by a cool spot, the RV amplitude is expected to be smaller in 
the IR regime. To test these two scenarios, we observed TW Hya 
with CRIRES, the CRyogenic high-resolution InfraRed Echelle 
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Fig. 4. CRIRES/H-band Radial-velocity points overplotted on 
the fitted Keplerian function shown in Table 1 . The error bars 
correspond to the quadratic addition of the errors bars on the RV 
standard and on TW Hya measurements. 



Spectrograph mounted on the VLT (iKaufl et al.l2006l) . To obtain 
the most accurate RV measurements, we selected a specific set- 
ting in the H band defined as 36/1/n (A c = 1.5884yum; AA~26/J.m) 
in the CRIRES manual. This setting is populated by deep and 
sharp CO2 lines that can be used as a simultaneous wavelength 
reference. The observations were done with a slit width of 0.2", 
using a typical ABBA nodding sequence, and a seeing > 0.8" 
in the optical - the objective was to reduce photocenter-induced 
RV errors while observing at the highest resolution available, es- 
timated as 100000. 

To derive the RV from the spec tra we used a cross- 
correlation method (Bar anne et alj|1996h . The spectra were cor- 
related simultaneou sly with a telluric m ask, built with the 
HITRAN database dRothman et alJll998l), and a s tellar mask, 
built from PHOENIX models dBarman et alj|2005l) . The atmo- 
spheric lines allow us to establish the zeropoint of the wave- 
length calibration. By subtracting the RV of the atmospheric 
lines from the target RV, we can correct from instrumental 
chan ges following th e same principle of the simultaneous ThAr, 
(e.g. lPepe et alj2007h . The whole procedure will be described in 
detail in a forthcoming paper (Figueira et aL, in prep.). 

TW Hya was observed five times between Feb 22 and 28 
2008. An RV standard (HD108309) was observed immediately 
after the target each night. We calculated the mean RV per night 
and assigned a lcr error bar corresponding to the standard de- 
viation of the RV in the nodding cycle. We verified that the 
measured scatter per night, between 20 a nd 50 m/s, matche d 
the expected uncertainties as estimated bv lBouchv et alj ((2001). 
However, the scatter of the 5 datapoints obtained on different 
nights is higher (~ 90 m/s), testifying the presence of systematic 
errors. Among others, all the observations in the first night were 
done with the guiding camera saturated. To remove these effects, 
and assuming that they affect both TW Hya and the RV standard 
in the same way, we calculated the relative velocity between TW 
Hya and the standard for each night. The rms of the five relative 
RV is 35 m/s, and 22 m/s if we only consider the 4 last points. 
The results are plotted in Fig. [4] and compared with the fitted 
orbit (Table [B- 



To estimate the match between the IR observations and the 
fitted orbit (Table Q]), we performed a Monte Carlo simulation. 
For each JD, we built an RV point by adding to the correspond- 
ing orbit value an observational error represented by a zero- 
centered Gaussian distribution with a fixed dispersion. We cal- 
culated the probability of obtaining the observed RV dispersion, 
hence the probability that our five points indeed belong to the fit- 
ted orbit. For the observed scatter of 35 m/s (the measured one) 
and 70 m/s (twice the observed one), the probability that our RV 
dispersion is drawn from the fitted orbit is always below 10~ 6 . 

Finally, we note that, as an important by-product of this 
work, we have obtained the most accurate RV measurements in 
the near-IR to date, showing the capabilities of CRIRES in this 
domain. 

3.3. Photometric observations 

TW Hya was monitored photometrically with the Euler tele- 
scope (La Silla, Chile) during 15 nights from 2008 February 2 
to March 13 under photometric conditions (except for March 8). 
Each night, a sequence of six exposures was observed at approxi- 
matively the same airmass to minimize any second-order extinc- 
tion effect. Observations were performed in the R-band and were 
defocused to ~ 5"to (1) minimize any systematics due to the 
intra/inter-pixels variations and (2) to allow exposure times of 
~ 50s in order to minimize the scintillation noise. The pointing 
was chosen to include TYC 7208-1066-1 and 3 other fainter iso- 
lated comparison stars in the EulerCAM field of view. Aperture 
photometry was used to extract the fluxes. Figure [5] shows the 
resulting lightcurve. 

TW Hya exhibited ~20% (peak to peak) flux variations 
during the run. The analysi s of the light curve using both the 
Lomb-Sc argle periodogram dScargld 19821) and the string-length 
method (Dworetskv l 19831) provides tentative periods of 6.1 and 
6.5 days, respectively. The origin of this period is not clear 
and the small time coverage does not allow us to draw fur- 
ther conclusions about this result. The comparison of this value 
with previous estimates does not help in interpreting our data, 
since the rotation period of TW Hya is not well-established 
and very different values (fro m 1.6 to 4.4 days) have been re- 
ported in the literature (e.g. [Mekkaden 1998; Batalha ^t al.l 
l2002ULawson & Crausdl2005HSetiawan et alj|2008l) . 

Interestingly, the photometric variations seem to be corre- 
lated with the surface of the CORALIE CCF (the CORALIE data 
and the photometric measurements were obtained on the same 
nights) - Fig|5j lower panel. Since the CORALIE CCF is a mea- 
surement of an average spectral line in the spectral region used, 
the measurement of its surface is rel ated to the veiling due to the 
hot spots on the stellar surface (e.g. lAndersen et alH l989). The 
correlation displayed in Fig. [5] would imply that the photomet- 
ric variations of TW Hya are dominated by hot spots, something 
that has already been rep orted in several studies (e.g. Mekkaden 
ll998fclBatalha et aLll2002h . The V-band amplitudes measured in 
these works range between 0.2-0.5 mag, that is, they are much 
larger than the 4% V-band photometric variations predicted by 
our cold-spot model. The fact that the photometric lightcurve 
cannot be fitted with the 3.56 days period derived from the spec- 
troscopic data (and assuming this as the true rotational period) 
can be related to the presence of short-lived hot spots (with typ- 
ical lifetimes of hours or days) t hat appear and dis appear hiding 
the true period of the source (see [Mekkaden 199j). 

We note that, although the effects of a hot spot on the 
RV variations are not tested here, we expect them to be sig- 
nificantly weaker than the ones induced by a cold spot. Note 
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Fig. 5. Upper panel: R-band lightcurve obtained for TWHya 
(filled symbols) and for the main reference star (open symbols). 
The fluxes are normalized. Lower panel: TW Hya photometric 
flux as a function of the surface of the CORALIE CCF. The best 
linear fit to the data is shown. 
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that their filling factors a re usually between 0.1% and 5% (e.g. 
iFernandez & Eiroalll996l) . that is, relatively small compared to 
the one of the dark spot from the model described in Sect. [2] 



4. Conclusions 

The analysis of new optical and infrared spectra of TW Hya al- 
lows us to conclude that the best explanation for the observed RV 
signal is the presence of a cold spot and not a hot Jupiter orbiting 
the star. First, we find a clear dependence of the optical RV am- 
plitude on the CCF mask used and, second, the IR RV curve is 
almost flat, with a small RV scatter (~35 m/s) that is inconsistent 
with the optical, orbital solution. In fact, we show that a simple 
cool spot model can reproduce the RV observations presented by 
SHL08. On the other hand, the photometric data are correlated 
with the surface of the CCF, which points towards a hot spot as 
the cause of the photometric variability, as reported in previous 
studies. 

Our result shows that searching for planets around young 
stars is highly demanding and cumbersome. Only the concerted 
use of the different techniques and instruments presented here 
were able to clarify the nature of the signal observed in TW Hya. 
It is, for instance, remarkable that no phase shifts were found in 
the RV signal between different epochs (Sect. 3.1), which shows 
that very stable cool spots exist around CTTSs. This must be 
taken with care when searching for planets around this kind of 
objects. 
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